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The Q dependence of the inelastic neutron scattering (INS) intensity of transitions within the 
ground-state spin multiplet of single-molecule magnets (SMMs) is considered. For these transitions, 
the Q dependence is related to the spin density map in the ground state, which in turn is governed 
by the Heisenberg exchange interactions in the cluster. This provides the possibility to infer the 
exchange-coupling constants from the Q dependence of the INS transitions within the spin ground 
state. The potential of this strategy is explored for the M = ±10 — > ±9 transition within the S = 10 
multiplet of the molecule Mni2 as an example. The Q dependence is calculated for powder as well 
as single-crystal Mni2 samples for various exchange-coupling situations discussed in the literature. 
The results are compared to literature data on a powder sample of Mni2 and to measurements 
on an oriented array of about 500 single- crystals of Mni2. The calculated Q dependence exhibits 
significant variation with the exchange-coupling constants, in particular for a single-crystal sample, 
but the experimental findings did not permit an unambiguous determination. However, although 
challenging, suitable experiments are within the reach of today's instruments. 

PACS numbers: 33.15.Kr, 71.70.-d, 75.10.Jm 



I. INTRODUCTION 

An important issue in the field of molecular nanomag- 
nets is the determination of the exchange-coupling con- 
stants between the metal centers in a cluster. In many of 
these molecules the exchange interactions are in fact the 
dominant terms in the spin Hamiltonian, and hence de- 
termine the nature of the ground state. The spin Hamil- 
tonian may be written as 



H 



^1 Jij^i 



Ha, 
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where the first term describes the Heisenberg ex- 
change interactions in the cluster (each exchange bond 
shall be counted only once) and the second term the 
weaker anisotropic contributions (single-ion anisotropy, 
anisotropic and antisymmetric exchange, dipole-dipole 
interactions). The first term will be also denoted as H ex . 

The single- molecule magnets (SMMs), such as the 
compound Mni2-acetate (or Mni2 in short), repre- 
sent important examples. The spectacular phenom- 
ena observed in these clusters, like the slow relax- 
ation of the magnetization or quantum tunneling of 
the magnetization^^ are intimately connected to a 
ground- state spin multiplet, which is characterized by a 
large spin S and a large anisotropy splitting of the easy- 
axis type. Obviously, the topography of the exchange- 
coupling constants is the key factor which controls 
the value of 5. However, it also has a strong impact 
on the magnitude of the anisotropy splitting in the spin 
ground state. On the one hand, it determines how ef- 
ficiently the microscopic anisotropy terms project onto 
the spin ground state, 5 and on the other hand, it gener- 
ates higher-order anisotropy terms via a mixing of spin 
multiplets (this mechanism is also called S mixing). ^iSaS 



Typically, the SMMs are modeled by the giant- 
spin Hamiltonian, which essentially describes the whole 
molecule by a single spin 5. The giant-spin Hamiltonian 
for Mni2, for instance, usually reads (in zero field) 



H S = D 



s 2 z -ls(s + i) 
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where S = 10, and D = -0.057 meV, B 40 -2.78 x 
10~ 6 meV, and B AA = -3.2 x 10~ 6 meV. 10 The dominat- 
ing D term splits the spin multiplet into ±M sublevels, 
with the M = ±5 level lowest in energy. It gives rise to a 
large energy barrier between the M = +S and M = —S 
states and hence the slow relaxation of the magnetiza- 
tion at low temperatures. The .840 and B44 terms are 
fourth-order anisotropy terms. 11 

The giant-spin Hamiltonian approach is phenomeno- 
logical in nature, but allows a very precise (in principle 
even exact) description of the energy eigenvalues of the 
ground- state multiplet, as well as their dependence on 
external parameters such as an applied magnetic field. 15 
One just has to choose the spin terms as appropriate (the 
possible spin terms are restricted by time reversal and 
cluster symmetry). The results of basically all experi- 
ments on SMMs which focus on the spin ground state, be 
it magnetization, specific heat, electron-spin resonance, 
nuclear magnetic resonance, Mossbauer spectroscopy, or 
else, can be very well described within this framework^ 

A notable exception is inelastic neutron scattering 
(INS). The INS selection rule AM = 0, ±1 permits tran- 
sitions between adjacent M levels, which results in a 
typical picket-fence INS spectrum with the M = ±5 — > 
±(5—1) transition at highest energy transfers j 16 i 17 While 
the giant-spin Hamiltonian works well for the transition 
energies, the Q dependence of the INS intensity cannot 
be described at all in this approach. In fact, the giant- 
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spin Hamiltonian approach predicts the INS intensity to 
be constant as function of Q ^ 16 ^ 17 while in reality it ex- 
hibits a typical oscillatory Q dependence due to the in- 
terference effects between the various spin centers in a 
cluster i 18 i 19 The Q dependence is hence a direct signa- 
ture of the many-spin nature of the ground-state spin 
multiplet. 7 The ability of INS to probe the many-spin 
character of wavefunctions is unrivaled by other tech- 
niques. 

This feature of INS suggests the following strategy 
to infer information on the exchange-coupling constants 
in a SMM. The INS intensity is governed by factors 
exp[iQ(R^ — Hj)](n\Si a \m) (m\Sj p\n) ^ i.e., depends on 
the geometrical arrangement of the spin centers in the 
cluster and the local spin transition matrix elements (Q 
is the momentum transfer vector, the position vec- 
tor of the ith spin center, \n) an eigenstate of the sys- 
tem, and a, (3 = x,y,z). The set of elements (n\Si a \m) 
may be called the spin transition map, 7 in analogy to 
the spin density map, (n\Si a \n). Via the Q vector in the 
exponential, the relative contribution of each of the lo- 
cal matrix elements (n\Si a \m) to the INS intensity can be 
modulated. This gives rise to interference factors and the 
characteristic Q dependence of the INS intensity. Hence, 
the Q dependence allows one to probe the topography of 
the spin transition map. 

Within the space of the ground-state spin multiplet 
of a SMM, however, the transition matrix elements 
(n\Si a \m) are proportional to the diagonal matrix ele- 
ments (n\Si a \n) in a first approximation, because of the 
Wigner-Eckart theorem: {tSM^^tSM') = C^ M ' x 
(tS\\Sj\\tS) = C^ M '/C^ M (rSM\S za \rSM), where 
(rS\\Si\\rS) is a reduced matrix element and C^f M 
are constants (essentially Wigner-3j symbols). In other 
words, the Q dependence of an INS transition within the 
ground-state multiplet is related to the spin density map, 
which in turn, as evident already from physical intuition, 
is determined by the exchange-coupling constants in the 
cluster. Hence, the Q dependence of the INS transitions 
between the ±M levels of the spin ground state of a SMM 
should allow one to retrieve the exchange-coupling con- 
stants. 

The novel concept in this approach should be noted. 
Usually, the coupling constants are determined by prob- 
ing, either via spectroscopic techniques or thermal excita- 
tion, the energies of the higher lying spin multiplets. The 
method suggested here, in contrast, retrieves the infor- 
mation exclusively from the ground-state multiplet. This 
could be, e.g., an advantage in cases where the higher ly- 
ing spin multiplets are not accessible to experiment for 
one or the other reason. 

This work aims at exploring this strategy with the ex- 
ample of Mni2. The Q dependence is calculated for vari- 
ous exchange-coupling topographies considered before in 
the literature, and compared to experiment. The two sit- 
uations of a powder and a single-crystal sample of Mni2 
are considered. In the powder case, the calculations are 



compared to the experimental data of Hennion et aZ.$2S 
for the single-crystal case data has been recorded on a 
quasi-single-crystal sample, which consisted of an ori- 
ented array of about 500 single crystals of Mni2. 

The outline of the manuscript is as follows. In the 
next section the relevant theoretical basics are presented. 
In Sec. Ill some exchange-coupling topographies are de- 
scribed and the reduced matrix elements or projection 
coefficients, respectively, calculated. In Sect. IV experi- 
mental details are given. Sections V and VI discuss the 
Q dependence of a powder and single-crystal SMM sam- 
ple, respectively. Section VII finally presents a conclu- 
sion. Some supporting information is provided in the 
Appendix. 



II. BASICS 

The differential neutron scattering cross section for 
spin clusters is given b y 21 i 22 



d 2 a _ m _e"^ 



dVtduo 
with 
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Here, C(Q,T) = {-fe 2 /m e c 2 )(k' /k) exp[-2W(Q, T)] (all 
symbols have the usual meaning), (3 = l/^^T), Z(T) is 
the partition function, Fi(Q) is the magnetic form fac- 
tor of the ith spin center, Q = k — k' is the transferred 
momentum, and R^- = R^ — Rj is the distance vector be- 
tween spin i and j. \n) denotes an eigenstate with energy 
E n of the microscopic spin Hamiltonian H. In this work, 
\n) is one of the 2S + 1 states of the ground-state mul- 
tiplet of a SMM. In the following, also the abbreviation 
la = Qa/Q is used. 

For a powder sample in zero magnetic field the INS 
intensity is obtained by averaging over all orientations of 
Q, which yields 19 

ij 

+32(QRij) E T f *( R ^') T f (SiSj)}, (5) 

q 

where jk is the spherical Bessel function of order 
Tq k \v) is the qth component of the spherical tensor of 
rank k constructed from the Cartesian vector v, and 
T^ 2 \SiSj) represents the tensor product [TW^) 
T^\Sj)] q 2 \ The ordered products Si a Sjp, which ap- 
pear in the explicit expression of Tf'iSiSj), stand for 
(n\Si a \m)(m\Sjp\ri). Equation ([5]) can be also written 
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more compactly ai 



,23 



(6) 



with the interference factors /^(Q, Rij) (fc = 0, 2 and 
|g| < fc) and the symmetrized spherical tensors Z7<^ 
(which are proportional to Re[T<^] for g > and 

Im[Tg^] for q < 0). Explicit expressions for Eqs. ([5j) 
and ([6]) are given in Refs. [7[ and [23[, respectively. 

Often it is desirable to describe the low-lying excita- 
tions of a spin cluster by an effective spin Hamiltonian 

H instead of the microscopic spin Hamiltonian H. The 
eigenstates of the effective spin Hamiltonian will be de- 
noted as |n), and the spin operators acting in the space of 

the effective spin Hamiltonian as Sjp (effective spin op- 
erators). Since H is supposed to yield the same energies 
as a distinction between E n and E n is not necessary. 
For a more detailed description of the relationships be- 
tween effective and microscopic spin Hamiltonian see e.g. 
Ref. Q. Using the effective spin Hamiltonian, the Q de- 
pendence of the INS intensity can be reproduced with 
a high degree of accuracy by a first-order perturbation 
theory approach.^^ 

In the present context of a SMM, where the excitations 
within the ground-state spin multiplet are of primary 
interest, the effective spin Hamiltonian H simply cor- 
responds to the well-known giant-spin Hamiltonian Hs- 
The effective spin operators Sjp are then related to the 
components of the total spin operator Sp via the projec- 
tion coefficients Ti(i), i.e., §i a = Fi(i)S a Z The matrix 
elements of the local spin operators §i a hence simply be- 
come 



(n\S ia \m) = ri(i)(n|5 a |ra). 



(7) 



The eigenstates \n) are obtained from the diagonalization 
of the giant-spin Hamiltoinan in the spin space \SM). 
The projection coefficients are determined as usual by 5 



ri(») 



(tS\\TW(S,)\\tS) 
(S\\TW(S)\\S) ' 



(8) 



where (rS\\T^ (Si)\\rS) is a reduced matrix element re- 
lated to the eigenstate \rSM) of the Heisenberg-exchange 
part H ex of the microscopic spin Hamiltonian, and 
(5||T( 1 )(5)||5> = ^(25 + 1)5(5 + 1) is the reduced ma- 
trix element in the spin space \SM). 

A subtle point shall be mentioned here. The described 
procedure essentially correspond to the strong-exchange 
limit, 5 i.e., the first-order perturbational treatment of the 
magnetic anisotropy. However, the strong-exchange limit 
does not yield all relevant spin terms in the giant-spin 
Hamiltonian correctly, such as for instance the fourth- 
order terms, because it neglects a mixing of the spin 
multiplets£i£&£ In principle, a similar statement is also 



true for the effective spin operators i However, it turns 
out that the mixing of the spin multiplets has a small ef- 
fect on the topography of the matrix elements (n\Si a \fh) 
even for rather large magnetic anisotropy, such that the 
Q dependence of the INS intensity is very well repro- 
duced in first order. This has been explicitly demon- 
strated in particular for the Mni 2 cluster, see Fig. 3 of 
Ref. [7| . Therefore one finds that the local spin operators 
are related to the total spin operator via factors (the pro- 
jection coefficients), which are determined completely by 
the Heisenberg interactions in the system. 

Insertion of Eq. ([7j) into the INS formula (|4|) , as appro- 
priate for a single-crystal material, yields 

Inm(Q) = F(Q) l a lp)Re (n\S a \fn)(fn\Sp\n) 

ap 

(9) 

with the interference factor 

.F(Q) =^F,(Q)F j (Q)cos(Q.R, j )ri(i)r 1 (j). (10) 

Because of the symmetry with respect to the interchange 
ij — > ji, the exponential is replace by the cosine and the 
real part, Re[], is introduced. It is interesting to note that 
the interference effects, described by ^(Q), factorize. 

For a powder sample, Eq. ([7j) should be inserted 
into, for instance, Eq. (|6]), which effectively replaces 



(ii) 



Uf'iSiSj) by r 1 (i)r 1 (j)/7^ ; (SS). The result is 
WQ) = ^^(Q)C/W(SS) 

kq 

with the interference factors 

f k9 (Q) = '£f%(Q,R ij )r 1 (i)r 1 (j). (12) 

In contrast to the single-crystal case, where the interfer- 
ence factor depends on the orientation of the transferred 
momentum vector, here the interference factors depend 
only on its magnitude. There is an instructive relation 
between Eqs. (|TQ|) and ([T2]) . which is outlined in the Ap- 
pendix. 

The giant-spin approach has been often used to ana- 
lyze INS data on SMMs i 16 i 17 However, as discussed in the 
introduction, it disregards the many-spin nature of the 
cluster; the interference factors are hence constants, and 
the predicted INS Q dependencies just flat. This fun- 
damental deficiency is overcome by Eqs. (|TQ|) and ([T2]h 
which reproduce Q dependencies accurately. They hence 
are much more preferable than the giant-spin approach. 
Importantly, for both the powder and single-crystal case 
the interference factors depend solely on the geometrical 
arrangement of the spin centers and the projection coeffi- 
cients on each site. They do not depend, for instance, on 
the particular form of the giant-spin Hamiltonian. Hence, 
as the structure of the cluster is known, the measured Q 
dependencies provide direct information on the projec- 
tion coefficients, and wherewith on the exchange-coupling 
constants. This is the basic idea of the present work, and 
Eqs. ([TO)) and ([T2j) are the explicit formulations of it. 
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FIG. 1: Schematic representation of the spin centers and 
exchange- coupling paths in Mni2. The spins 1-4 are denoted 
as core and the spins 5-12 as the crown. The right graphic 
shows the molecular structure of Mni2 viewed along the crys- 
tallographic c axis (only the Mn and the bridging oxygen 
atoms are shown). The orientations of the crystallographic 
a and b axes of the tetragonal unit cell are indicated. The 
core and crown are rotated with respect to the a and b axes 
such that the axis through spins 1 and 3 is rotated by 12.33° 
and that through spins 7 and 11 by 10.91°. 



III. PROJECTION COEFFICIENTS 

As shown in the previous section, the evaluation of the 
INS Q dependence involves the geometrical structure of 
the cluster and the projection coefficients. The molecu- 
lar structure of Mni2 is precisely know from x-ray crys- 
tallography (Mni2 crystallizes in space group I4 and ex- 
hibits a molecular S4 symmetry axis, see also Fig. [Q)^ i n 
this section, the projection coefficients for Mni2 will be 
determined for various exchange-coupling topographies 
considered before in the literature. 

The exchange-coupling graph appropriate for Mni2 is 
presented in Fig. [TJ It involves the four exchange con- 
stants Ji, J2, J3, and J4. Due to the S4 symmetry of the 
cluster, the projection coefficients assume three different 
values, namely 

Ti(A) = imi) = ri(2) = ri(3) = ri(4), (13) 
r-i(B) = r 1 (6) = r 1 (8) = r 1 (io) = r 1 (i2), (14) 
rye) = ri(5) = r-i(7) = r\(9) = r-i(ii). (15) 

The indices A, £?, and C will be used also to 
notate the three sublattices consisting of the spins 
{1,2,3,4}, {6,8,10,12}, and {5,7,9,11}, respectively. 
Since T 1 (A) + T 1 (B) + Ti(C) = 1/4 the projection co- 
efficients are not independent of each other, and the Q 
dependence of the INS intensity is effectively governed 
by only two independent parameters. 

The first idea on the coupling situation in Mni2 came 
from magneto-chemical considerations, which suggested 
that J\ is antiferromagnetic and much larger than the 
other three coupling constants. 25 Hence, the spins of e.g. 
the Mn ions 1 and 7 couple to an overall spin of 1/2, 
and similar for the spins 2&9, 3&11, and 4&5. In order 
to arrive at aS = 10 ground state, these spins are then 
coupled ferromagnetically to the remaining four Mn(III) 



spins 6, 8, 10, and 12. Following Ref. [26], this scheme 
will be called the "Florentine" coupling scheme. Since 
the ground-state spin wavefunction is constructed here by 
a successive coupling of spins, the projection coefficients 
can be calculated easily using the Wigner-Eckart theorem 
and Racah formalism. 5 

A second scheme consists of coupling the eight Mn(III) 
spins ferromagnetically to yield an overall spin Sbc = 
16, to couple the four Mn(IV) spins ferromagnetically 
to yield a spin Sa = 6, and then to couple these two 
spins antiferromagnetically to yield S = Sbc — Sa — 
2Q^26 rpfiis coupling scheme is based on the rational that 
J\ and J 2 are both antiferromagnetic, while J3 and J4 
are assumed to be ferromagnetic. It will be called the 
"A+(B+C)" scheme. The projection coefficients again 
can be calculated easily using the Wigner-Eckart theorem 
and Racah formalism. 

A third scheme consists of coupling the four spins of 
each sublattice A, £?, C ferromagnetically, then to cou- 
ple the spins of sublattices A and B antiferromagneti- 
cally to yield Sab = 4, and to finally couple this spin 
with that of sublattice C antiferromagnetically to yield 
S = Sab — Sc = 10. This coupling scheme is based on 
the rational that J\ and J2 are both antiferromagnetic 
and dominating, such that the frustration induced by J3 
and J 4 is negligible. It will be called "(A+B)+C". The 
calculation of the projection coefficients proceeds as for 
the previous two schemes. 

The above three schemes approximate the wavefunc- 
tion of the ground-state spin multiplet by a spin wave- 
function \vSM), where v denotes intermediate spin quan- 
tum numbers. The associated spin-coupling schemes 
were suggested by the different assumptions concerning 
the exchange-coupling topography. A more penetrating 
approach is of course to specify the values of the coupling 
constants J\ , . . . , J4 and to calculate the ground-state 
wavefunction numerically by exact diagonalization. Such 
an approach has been pursued for Mni2 by three groups. 
Raghu et al. suggested the set J\ = — 215 K, J 2 = J3 = 
—85 K, J4 = 64.5 K*2£ The projection coefficients can be 
determined from the spin densities given in their work. 
Regnault et al. deduced the values J\ = —119 K, J 2 — 
-118 K, J 3 = 8 K, J 4 = -23 K. 28 No information is given 
which would allow a determination of the projection co- 
efficients. Hence we calculated them numerically using 
a sparse-matrix iterative subspace techniquei2£ Finally, 
Chaboussant et al. deduced the values J\ — —67.2 K, J2 
= -61.8 K, J 3 = -7.8 K, J 4 = -5.6 K. 30 ^ 31 Here too the 
projection coefficients were calculated by us numerically. 
The exchange-coupling constants were also evaluated by 
LDA+U and DFT calculations ,22*22 but these additional 
sets do not change the conclusions of this work and are 
hence not considered. The works Refs. [34,35,36], which 
start from an effective 8-spin model for Mni2, are men- 
tioned for completeness. 

The projection coefficients of the six exchange-coupling 
topographies just described are compiled in Table H For 
all cases, T\(A) is negative while T\(B) and T\(C) are 
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TABLE I: Projection coefficients for the exchange-coupling 
topographies discussed in the text. 



scheme 


ri(A) 




IMC) 


Florentine 


-0.05 


0.2 


0.1 


A+(B+C) 


-0.13636 


0.19318 


0.19318 


(A+B)+C 


-0.1 


0.15 


0.2 


Raghu et al. 


-0.09019 


0.18103 


0.15916 


Regnault et al. 


-0.12472 


0.18376 


0.19096 


Chaboussant et al. 


-0.12121 


0.18398 


0.18723 



positive. This reflects the antiferromagnetic alignment of 
the spins of the Mn(IV) core with respect to the spins on 
the Mn(III) crown. The different exchange-coupling to- 
pographies clearly express themselves in different projec- 
tion coefficients. However, the variation is not very pro- 
nounced, with the exception of the Florentine coupling 
scheme. For this scheme T\{A) and T\{C) are signifi- 
cantly smaller in magnitude, which reflects the assump- 
tion of a dominant J\ in this scheme. 



IV. EXPERIMENTS 

In the next sections, the Q dependence of the INS tran- 
sition between the M = ±10 and M = ±9 levels of the 
S = 10 ground-state multiplet of Mni2 will be consid- 
ered. It occurs at an energy transfer of 1.24 meV, and is 
the strongest inelastic peak at liquid- He temperatures. 

The best measurement so far of the Q dependence of 
this transition in powder samples of Mni2 has been re- 
ported by Hennion et al.£& These authors studied a par- 
tially deuterated Mni2 sample on the triple-axis spec- 
trometers 4F1 and IT at the Laboratoire Leon Brillouin, 
Saclay, France, at a temperature of 1.55 K. The calcula- 
tion of the Q dependence for powder Mni2, which will be 
presented in the next section, hence will be compared to 
this data set. 

In order to determine the Q dependence of a single- 
crystal sample, we performed measurements on an ori- 
ented array of about 500 non-deuterated single crystals of 
Mni2. The crystals were synthesized following literature 
procedures. 24 The dimensions of the needle-shaped crys- 
tals were approximately 5x0.5x0.5 mm 3 , and the masses 
ca. 2 mg. Mni2 crystallizes in the space group 14, hence 
the needle axis coincides with the magnetic anisotropy 
axis z. The oriented array of single crystals was obtained 
by placing crystals in long, narrow grooves milled into 
aluminium platelets, such that the needle (= z) axes were 
aligned along the grooves. 37 A total of 15 such platelets 
which each held about 30—35 single crystals were stacked 
in an aluminium container and sealed therein. The INS 
intensity was measured on the time-of-flight spectrome- 
ter IN5 at the Institut Laue-Langevin, Grenoble, France. 
The initial neutron wavelength was set to 5.9 A; the res- 
olution at the elastic line was 60 /ieV. The sample was 



inserted into an orange cryostat, permitting a sample 
temperature of 1.5 K. The data were corrected for the 
contribution of the sample container and the detector ef- 
ficiency using a vanadium standard. 

The alignment of the z axis of the quasi-single-crystal 
sample with respect to the transferred momentum vec- 
tor Q crucially influences the measured Q dependence. 
In simulations, details of the experiment have hence to 
be considered. The laboratory frame XY Z is chosen such 
that the X axis points along the wave vector of the inci- 
dent horizontal neutrons, ex = k/fc, and the Z axis into 
the sky. In the experimental configuration of the IN5 
spectrometer, the Y axis is then directed towards the 
side of the detector banks. The Q vector lies in the XY 
scattering plane and is described by its magnitude Q and 
angle ipQ, Q = Q(cos (pQ, sin <^q, 0), with Qy < 0. The 
calculation of matrix elements, however, is most conve- 
niently done in the xyz frame of the molecular magnetic 
axes. As will be shown in the next section, the cluster 
can be treated as uniaxial, and the orientation of the 
sample hence described by the polar angles 6 and cp of its 
z axis in the XYZ frame. This procedure is not exactly 
correct, because the cluster coordinates are not invariant 
under rotations around z. We found, however, that this 
subtlety does not affect the conclusions of this work and 
hence do not exploit it further here. In our calculations 
the x, y axes were along the crystallographic a, b direc- 
tions, see Fig. [U The z axis in the XYZ frame is gener- 
ated by two rotations, z = Rz((p)Ry(0)Z. The transfor- 
mation of a vector from the XYZ to the xyz frame is thus 
Ry 1 {0)R~ 1 {Lp), and for Q in the xyz frame one obtains 
Q = Q(cos#cos((/?q —<p), sin((/?Q - <p), sin# cos((p Q -cp)). 
The experimental data was taken for four different config- 
urations, namely with the z axis along (0,0,1), (1,0,0), 
(0, 1, 0), and (—1,1,0), where (X, Y, Z) denotes a vector 
in the XY Z frame. 

For comparison of experiment and calculation, it also 
has to be considered that the IN5 spectrometer measures 
the intensity as function of the direction k'/fc', and not Q. 
The data analysis software LAMP, 39 however, which was 
used for data reduction, yields the intensity as function 
of Q (which is related to cpq via the scattering triangle). 
Hence, the experimental Q dependencies include an im- 
plicit variation of (pQ, i.e., the curves do not correspond 
to Q scans as they would be obtained on, e.g., triple-axis 
instruments. 



V. Q DEPENDENCE OF A POWDER SMM 
SAMPLE 

In this section, the Q dependence of the M = ±10 — > 
±9 transition of a powder sample of Mni2 will be cal- 
culated for the six exchange topographies described in 
Sec. Ill, and compared to the experimental results of 
Hennion et al. 20 In order to proceed, the matrix ele- 
ments (n\S a \fh) have to be evaluated. Here \n) denotes 
one of the two states of the M = ±10 level, and \m) 
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one of the two states of the excited M = ±9 level (the 
difference between \n) and \n) is recalled, Sec. II). 

Some points need clarification. First, a zero-field situ- 
ation is considered. Hence, the wavefunctions may be 
chosen as real. The M = ±10 and M = ±9 lev- 
els each consist of two states, which shall be denoted 
as |0)± and |1)±, respectively. Since for these levels 
the tunnel splitting is way to small to be observed in 
INS, the intensity of the M = ±10 — > ±9 transition 
in fact originates from four transitions. This is eas- 
ily accounted for in the equations by considering each 
square of matrix elements to correspond to the sum 
Z) M i/=±(°U'S'a|l)i/(l|i/S , ^|0) M . Second, the eigenstates of 
the giant-spin Hamiltonian are not eigenstates of S z be- 
cause of the presence of terms such as B^O\. 11 They 
hence cannot be labeled by M (we nevertheless use the 
notion of, e.g., a M = ±10 level, its correct meaning 
should be obvious). The four lowest-lying states in fact 
are very well described by |0)± = (|10) ± | - 10))/>/2 and 
|1)± = (|9) ± | - 9))/V2. This means (i) that the matrix 
elements for the operator S z are essentially zero, and 
(ii) that the Mni2 molecule can be considered as mag- 
netically uniaxial (in fact S MI/=± (0| At 5 QJ |l) I/ (l| I/ S^|0) At = 

£m=±io,M'=±9 ( M I Sa I M') (M f | Sp | M) ) . These consider- 
ations are generally valid for the states of a SMM at the 
bottom of the energy barrier (in zero magnetic field). 

With the choice of real wavefunctions, (n\S z \fh) = 0, 
and uniaxiality, only Z7q°^ and in Eq. ([6]) are non- 
zero with Uq 2 ^ = —I/^/EUq ^. Hence one obtains 



Inm{Q) = \U { ] F P {Q) 



(16) 



with the interference factor 



^p(Q) = 5^F»(Q)Fi(Q)ri(i)ri0-) 



'\jo(QRi; 



3R 2 

h(QHu) l:hZ 



E>2 

-£].(17) 



This is a rather general result. It holds for the INS tran- 
sitions in a powdered SMM at the bottom of the energy 
barrier [it also implies that the Q dependence does not 
depend on M, only for states near the top of the barrier 
it may display deviations from !F P {Q) due to the effects 
of S z ]. Interestingly, the interference factor again factor- 
izes. The single-spin approach yields I nm (Q) oc Uq°\ i.e., 
Fp(Q) = const iLiSaii Comparison with Eq. (fTB|) hence 
again nicely demonstrates the effect of the many-spin na- 
ture of a real spin cluster, which is expressed here by the 
interference factor J- p (Q). 

The Q dependencies calculated from Eq. (jTTJ) are 
shown in Fig. Efa) for the projection coefficients given 
in Table I. The form factors Fi{Q) were estimated with 
the standard analytical approximations. 38 The non-zero 
value for Q — > is typical for a AS = transition. 22 
After a steep drop with a minimum at about 0.5 A -1 , 
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Raghu et al. 
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FIG. 2: (Color online) (a) Calculated Q dependencies of the 
intensity of the M — ±10 — > ±9 transition in a powder sample 
of Mni2 for the six exchange- coupling topographies discussed 
in the text. Curves are labeled in the order of the height 
of the maximum at 1 A -1 , (b) Experimental data of Hen- 
nion et al (reproduced from Ref. |20I|V For comparison also 
the calculated curves for the Chaboussant et al and Floren- 
tine coupling schemes are reproduced (the Chaboussant et al 
curve was scaled such as to reproduce the data, the Florentine 
curve was then scaled by the same factor). 



the curves exhibit a first pronounced maximum at about 
1.0 A -1 , display a broad feature between ca. 1.5 and 
2.8 A -1 , followed by further oscillations at higher Q val- 
ues. Apparently, although significant variation is seen, 
the Q dependence of powder samples does not depend 
very sensitively on the exchange-coupling topography. 
The most significant effect is displayed by the height of 
the first maximum, which is most reasonably measured 
with respect to the two minima to the left and right. A 
qualitative comparison of Fig. Efa) with Table I reveals 
that the height of the first maximum is correlated to the 
magnitude of Ti(A). However, presumably it would re- 
quire very precise experiments to measure this height re- 
liably. The broad feature in the range 1.5-2.8 A" 1 looks 
very similar for all considered cases, with the exception of 
the Florentine coupling scheme, for which it exhibits two 
rather well resolved maxima at about 1.8 and 2.4 A -1 . 

Figure EJb) shows the experimental data from Ref. [50 
for comparison. All calculated Q dependencies reproduce 
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the experimental curve more or less (the Chaboussant et 
al. curve has been also drawn as a representative for the 
five schemes other than the Florentine scheme, which all 
give similar agreement). The Florentine coupling scheme, 
however, shows the least agreement, see Fig. [2fb). First, 
the experimental height of the first maximum is clearly 
larger than calculated (in particular if one considers that 
e.g. incoherent scattering due to hydrogens should re- 
duce the experimentally observed height). Second, the 
broad feature seems to be better reproduced by the other 
curves. 

Hence, as a summary, the powder Q dependence shows 
some variation for different exchange-coupling topogra- 
phies, but they are not easily differentiated in the exper- 
iment. However, it is fair to conclude that the Florentine 
coupling scheme is not consistent with the data, and can 
be ruled out (as has been found also befor o 27 i 28 i 30 i 31 ). 
Hence, the powder Q dependence has some potential for 
differentiating between very different coupling topogra- 
phies (it should be recalled that the Florentine scheme is 
distinguished from the others by its very dominant J\). 



VI. Q DEPENDENCE OF A SINGLE-CRYSTAL 
SMM SAMPLE 




0.5 1.0 1.5 



Q (A 1 ) 



0.5 1.0 1.5 

Q (A 1 ) 



FIG. 3: Measured Q dependencies of the intensity of the M — 
±10 — > ±9 transition in a quasi-single-crystal sample of Mni2 
for the four orientations of the z axis indicated in each panel. 
The normalization is arbitrary but identical for each curve; 
the curves hence correctly reflect relative intensities. 



In this section, the Q dependence of the M = ±10 — > 
±9 transition of a single-crystal sample of Mni2 will be 
considered. In this case, the INS intensity does not only 
depend on the magnitude of Q, as in the powder case, but 
also on its orientation with respect to the sample. The 
orientation of Q enters both Eqs. ([9]) and ([TO]) . Equa- 
tion (j9]), however, can be further simplified. With real 
wavefunctions, the matrix elements of S x and S z are real, 
while those of S y are imaginary. Hence the sum over a, 
(3 reduces to 

WQ) = ^(Q)E(l-^)|(n|5 a |m)| 2 -2U, 

a 

x(n\S x \fh)(n\S z \m)]. (18) 

Further considering that the matrix elements of S z are 
basically zero and that the cluster can be treated as 
magnetically uniaxial, as discussed in section V, one ar- 
rives at the approximation i nm (Q) = J r (Q)(2 — l\ — 
ll)\(n\S x \m)\\ 

Figure [3] shows the experimental results for the inten- 
sity of the M = ±10 — » ±9 transition in our quasi-single- 
crystal sample of Mni2 for the four orientations with z 
parallel to (1,0,0), (-1,1,0), (0,1,0) and (0,0,1), re- 
spectively (the vectors refer to the XY Z frame). Clearly, 
the Q dependence of the INS intensity depends very 
markedly on the orientation of the Q vector with respect 
to the sample. The calculated Q dependencies are pre- 
sented in Figured! for the six exchange-coupling topogra- 
phies listed in Table I. 

Comparison of Figs. [3] and [4] reveals that the experi- 
mental and theoretical curves exhibit very similar shapes 



for each orientation. In particular, the positions of the 
minima and maxima are in good agreement. For the 
(0,0,1) direction, the experimental data shows a maxi- 
mum at about 0.95 A -1 , which is embraced by two min- 
ima at about 0.5 and 1.4 A -1 with a slight upturn to- 
wards higher Q. In the (1,0,0) direction, the experi- 
ment exhibits a steep drop towards a minimum at about 
0.65 A -1 , then a broad maximum at about 1.2 A -1 
which is followed by a slow drop at higher Q values. The 
(—1,1,0) curve is characterized by a broad asymmetric 
feature with a maximum at about 0.8 A -1 and a steeper 
flank at the side of higher Q values. In the (0, 1,0) di- 
rection, finally, the experiment shows a peak at about 
1.1 A -1 with indications of two minima at about 0.6 and 
1.5 A -1 . The theoretical curves reproduce these features 
very well. 

However, also differences can be noted. In particu- 
lar, the scattering intensities at the minima are signifi- 
cantly higher than predicted by theory. In fact, many 
of the minima are predicted to have almost zero inten- 
sity. This is not unexpected, it is a rather typical obser- 
vation in non-deuterated samples and due to the large 
incoherent scattering of the hydrogen atoms i 23 i 4Q i 41 Fur- 
thermore, while the relative intensities for the three direc- 
tions (1,0,0), (—1,1,0), and (0,1,0) are in accord with 
the theoretical trend, the (0,0, 1) curve is almost a fac- 
tor of 2 more intense than predicted. This could be an 
effect of the sample geometry, which is asymmetric in 
the sense that a neutron beam along the z axis has to 
traverse about 1.5 cm of Mni2 while a neutron beam in 
the xy plane traverses at most about 0.8 cm. This sug- 
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FIG. 4: (Color online) Calculated Q dependencies of the in- 
tensity of the M = ±10 — > ±9 transition in a single crystal of 
Mni2 for the orientations of the z axis indicated in each panel. 
Each panel shows the results for the six exchange- coupling to- 
pographies listed in Table I. 



gests that absorption is weaker in the (0,0, 1) direction, 
consistent with the observation. 

The theoretical curves in the (—1,1,0) and (0, 1, 0) di- 
rections are basically insensitive to the exchange-coupling 
topography, while for the height of the maximum in the 
(0,0,1) and (1,0,0) directions the same trend as found 
for the powder Q dependence, Fig. [2j is observed. Based 
on the very weak intensity expected for the Florentine 
coupling scheme, which is in contrast with the experi- 
ment, this scheme again is ruled out. 

Besides the above conclusions we feel that the accuracy 
of the present experiment is not such as to allow further 
conclusions. At this point it is perfectly obvious to ask: 
how well can the projection coefficients be retrieved from 
the Q dependence and what would be an appropriate 
procedure. In order to explore this, the calculated INS 
intensity for the various exchange-coupling topographies 
for Mni2 is shown in Fig. [5] for the (0,0, 1) direction as 
function of a momentum transfer in the XY plane [the 
(0,0, 1) curves in the previous Fig. [4] correspond to cuts 
along the circle (Q x - 1.065) 2 + Q 2 y = (0.72) 2 ]. Evi- 
dently, the single-crystal INS intensity exhibits features 
of significant intensity at certain special points, which 
are related to the molecular structure (vide infra). The 
intensity of these "towers", however, varies significantly 
with the exchange-coupling topography. A clear trend 
of the intensities is seen along the sequence A+(B+C), 
Regnault et a/., Chaboussant et a/., (A+B)+C, Raghu 
et al., and Florentine, which is the same trend noted be- 
fore for the height of the maximum at about 1 A -1 [see 
Figs.E^a) and [J. 



Regnault et al. 
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4 -3 -2-10 1 
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FIG. 5: (Color online) Calculated INS intensity of the M = 
±10 — > ±9 transition in a single crystal of Mni2 for the (0,0,1) 
orientation as function of the momentum transfer in x and y 
direction (Q z = 0). The INS intensity is represented by the 
color [dark gray (blue online) = 0, light gray (red online) = 
maximum]. Each panel shows the result for the indicated 
exchange- coupling topography. 



A second, more striking trend can be observed. At 
Q w 1 A -1 there is a circle of intensity, which exhibits 
four minima and maxima (if one goes around the full 
360°). While for the (A+B)+C scheme the left maxi- 
mum (in the quadrant Q y < 0, Q x > 0) is at about 
—27° from the y axis, it rotates to about +7° in the 
sequence of (A+B)+C, Regnault et a/., Chaboussant et 
al, A+(B+C), Raghu et al, and Florentine [the Floren- 
tine scheme in fact has its maxima where the (A+B)+C 
scheme has its minima, and vice versa]. Observation of 
the locations of the minima and maxima on the 1 A -1 
circle would be hence very interesting. 

A qualitative understanding of the intensity pattern 
can be obtained if one disregards the Q dependence of 
the form factors Fi(Q), i.e., assumes Fi(Q) = 1. The 
cos(Q • Rij) factors in J-(Q), Eq. (fT0]h can be identically 
replaced by e zQ RiJ [the same sum but with sin(Q • R^ ) 
equals zero]. Then one finds 

f(Q) = J dRe iQ - R S(R), (19) 

g(R) = ^ri(i)riOXR-Rii), (20) 

i.e., T(Q) is just the Fourier transform of G(H). It should 
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FIG. 6: (Color online) ^(Q), Eq. (HUJ , of Mm 2 in the (0,0,1) 
orientation as function of the momentum transfer in x and 
y direction (Q z = 0). The value of ^(Q) is represented by 
the color [dark gray (blue online) = 0, light gray (red online) 
= maximum]. The panels show the results for the indicated 
situations discussed in the text. 



be noted that the INS intensity 7 nm (Q), Eq. ([9]), in gen- 
eral cannot be Fourier transformed because of the Q de- 
pendence of the sum. Exceptions are special cases such 
as uniaxial clusters with Q in the xy plane. Then the 
sum is a constant, see also Eq. ([T8]) . and also 7 nm (Q) is 
Fourier transformable [if Fi(Q) = 1]. Since *F(Q) is the 
Fourier transform of the distance vectors R^/ (and not 
of the position vectors as in diffraction), it exhibits 
maxima at Q • ~ 27rn, which are, however, severely 
broadened since only a small number of metal centers is 
involved (n is an integer here). 

It is revealing to plot ^(Q) for the cases 
{ri(A),ri(B),ri(C)} = {1,1,1} ("crown+core"), 
{-1,1,1} ("crown-core"), {0,1,1} ("crown only"), and 
{1,0,0} ("core only"), as presented in Fig. [6j The 
crown+core case represents the direct Fourier transfor- 
mation of the "density of distances" . 5(H — R^-) and 
mimics a fully ferromagnetic situation. The crown-core 
case in contrast mimics an antiferromagnetic alignment 
of the core spins with respect to the crown spins. The 
other two cases focus on the core and crown spins, re- 
spectively, i.e., neglect the interferences between core and 
crown. 

The core-only plot demonstrates that indeed broad- 
ened peaks appear at Q • R^ « 2i\n in ^(Q). The com- 
parison of the crown+core and crown-core plots shows 
that they are rather complementary in the sense that at 
the special Q points the former has maximum intensity 
where the latter has minimum intensity, and vice versa. 
Most interestingly, the intensity circle at Q « 1 A -1 is 
entirely absent in the crown+core plot while it is rather 
strong in the crown-core plot. It hence directly reflects 
the antiferromagnetic spin alignment of the core and 



crown spins. The crown-only plot shows the contribu- 
tion from the interference on the crown. It contributes 
only very weakly to the intensity at Q ~ 1 A -1 . 

This discussion shows that each feature visible in one 
of the plots of Fig. [5] can be related to certain pairs 
of spin centers, and its amplitude to Ti(i)Ti(j). One ap- 
propriate strategy to determine the values of Ti(i)Ti(j) 
would be thus to measure the INS intensity at the Q 
points where Q • R^- = 27rn. Since most features ap- 
pear for Q > 2 A" 1 , it is now also obvious that the 
experiment presented in this section, which was limited 
to 0.4 < Q < 1.6 A -1 , could not be very successful in 
discriminating the different exchange-coupling topogra- 
phies. 



VII. CONCLUSIONS 

In this work the possibility of extracting the exchange- 
coupling constants of SMMs from the Q dependence of 
the INS intensity of transitions within the ground-state 
spin multiplet has been explored, with Mni2 as an ex- 
ample. The calculated Q dependence of powder SMM 
samples was found to show significant variations for dif- 
ferent exchange-coupling situations, but in order to dis- 
tinguish these experimentally, very precise measurements 
are required. The powder Q dependence has a potential 
to discriminate between very different exchange-coupling 
situations. In the case of Mni2, the Florentine coupling 
scheme, which assumes a very dominating J\ interaction, 
could be discarded. The Q dependence of single-crystal 
SMM samples has a higher potential, because the addi- 
tional dependence of the INS intensity on the orientation 
of the momentum transfer vector Q allows for more in- 
dependent measurements. The presented experimental 
results on a quasi-single-crystal Mni2 sample, which con- 
sisted of an oriented array of about 500 single crystals, 
however, did not allow us to infer information beyond 
that already concluded from the powder Q dependence. 
On the one hand, this is due to the obvious limitations 
of the experiments. On the other hand, a qualitative dis- 
cussion of the situation shows that these experiments did 
not probe all the relevant Q space and indicates a better 
measurement procedure. 

As a summary, deducing the exchange-coupling con- 
stants in SMMs from the INS Q dependence seems to be 
a viable strategy, but is experimentally challenging. It 
requires precise work in order to arrive at quantitative 
conclusions concerning the coupling constants. For in- 
stance, highly deuterated samples seem to be mandatory. 
Also, absorption corrections due to the actual sample ge- 
ometry probably have to be employed, very much as in 
X-ray crystallography. However, although these issues 
clearly need further development, they are within reach 
of today's spectrometers. 

The key point in the approach is the relationship be- 
tween the exchange-coupling constants and the spin den- 
sity map. For the SMMs it was noted in this work that 
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the spin density map can be determined from the Q de- 
pendence of the INS transitions within the spin ground 
state of the SMM. However, a similar approach is clearly 
also possible with other experimental techniques, which 
allow one to measure the spin density map, such as elastic 
neutron scattering or nuclear magnetic resonance. 
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APPENDIX A 

In this Appendix, a relationship between the INS 
formulas for single-crystal and powder SMMs, i.e., 
Eqs. (2D, (TO) and (fTT|) . (TH2j) . respectively, shall be dis- 
cussed. As mentioned already in Sec. VI, the interference 
factor ^(Q) of a single-crystal SMM, Eq. ([TO]) , can be 
written identically with cos(Q • Hij) replaced by e zQ RiJ . 
Expanding the exponential in spherical harmonics, J~(Q) 
can be rewritten as 

.F(Q) = 47r^(i) fc (-l)«y«(Q)C fc ^(Q), (Al) 

kq 

where 

r kq {Q) = ^/^R^r^r^-), (A2) 

/ fc «(Q,Ry) = F i (Q)F i (Q)j fc (Qi? ij )T ( | fc )(R ii ). (A3) 

Here, k may be restricted to even values, k = 0, 2, 4, — 
Yq ^ (V) denotes the qth component of the spherical har- 
monics of rank k related to the Cartesian vector V. It is 
proportional to the irreducible spherical tensor Tq k ^ (V) , 

Y q ik) (V) = C k T^\v), which also defines the propor- 
tionality constant Ck q - 

This is the main finding of this Appendix: Equa- 
tions (|A2p and (|A3ft are identical to the results for the 
interference factors of powder SMM samples — if one ex- 
pands I n m(Q) m terms of the spherical tensors Tq^ and 
not the symmetrized spherical tensors as in Eq. (fTTj). 

This point deserves comment. Equations (|A3ft and ([T2]) 
look exactly identical, but in fact differ in the definition of 
the interference factors f kq (Q, Rij): Equation ([T2]) is de- 
rived from Eq. ([6]) , which in turn is obtained by rewriting 
the general INS formula for powder samples, Eq. (5]), in 

(k) 

terms of the symmetrized spherical tensors Uq . Hence, 
the expressions for the f kq (Q,Hij) contain the tensors 

Uq k \llij). 7 An expansion in terms of Uq^ is most con- 
venient for explicit calculations. The general powder for- 
mula Eq. ([5]), however, also can be expressed in terms of 



the spherical tensors Tq ' . The resulting equations look 
exactly the same as that for Uq k \ but with all Uq^ re- 
placed by Tq k>} [for instance, replacing Tq k>} by Uq^ in 
Eq. (A3l) directly yields the equations (B2) of Ref. 0]. 

(k) 

In this Appendix, the Tq -expansion is preferred because 
the algebra is then more obvious, but in the expressions 
all Tq k ^ can be replaced identically by Uq k \ In this sense, 
Eqs. (|A2p and (A3]) are identical to the results for powder 
SMM samples, such that the above equations establish 
the sought-after relationship. 

The relation Eq. (|A1|) is convenient. It separates the 
dependence of the INS intensity on the orientation of 
the Q vector from the dependence on the magnitude 
Q and the cluster properties [which enter via R^ and 
ri(i)]. Furthermore, it establishes the connection be- 
tween the single-crystal and powder interference factors 
(for the powder INS intensity, however, only the factors 
with k = 0, 2 are relevant, while in the single-crystal case 
no such restriction exists). 

Using the orthogonality relation for spherical harmon- 
ics, one obtains 

J ^^(Q)FW(Q) = (i) k C k F k *{Q). (A4) 

Hence, integrating the interference factor ^(Q) over all 
orientations of Q for a given magnitude Q allows one to 
extract the interference factors F kq (Q). 

If one assumes Fi(Q) = 1 for the moment, and uses the 
closure relation of spherical Bessel functions, one obtains 

J dQQ 2 T^\Q)j k {QR) = 
^ ^r 1 (i)r 1 (j)T«(R ij )5(i? - Rij), (A5) 

ij 

where R is an arbitrary number [the result also holds 
for the weaker condition of Fi{Q) = F(Q) for all i if 
J rkq (Q)/F(Q) is considered instead of F kq {Q)\. Hence, 
integrating JF^Q) over Q using a value of R which 
matches one of the metal distances Rij allows one to 
project out the value of Ti(i)Ti(j) for the pair [to 
be precise, one has to assume that if R^ is equivalent 
for several pairs (ij) that then also Ti(i)Ti(j) is equiva- 
lent for these pairs, which is reasonable]. If one combines 

Eqs. (AH) and (A5]) and sums over Y, kq T q k) *(R>ij), then 
one recovers the Fourier transform of ^(Q), Eq. ([20]) . 

A final comment: Obviously, very similar mathematics 
can be performed in order to establish relations between 
the INS formulas Eq. (H) and Eq. (6J, which are valid for 
general spin clusters. Following the above lines, the cor- 
responding equations are simple to obtain. They hence 
shall not be reproduced here. 
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